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ABSTRACT 

This work presents a detailed analysis of the overall flow structure and unique features of the inner 
region of the tilted disk simulations described in Fragile et al. (2007). The primary new feature 
identified in the main disk body is a latitude-dependent radial epicyclic motion driven by pressure 
gradients attributable to the gravitomagnetic warping of the disk. The induced motion of the gas is 
coherent over the scale of the entire disk and is fast enough that it could be observable in features 
such as relativistic iron lines. The eccentricity of the associated fluid element trajectories increases 
with decreasing radius, leading to a crowding of orbit trajectories near their apocenters. This results 
in a local density enhancement akin to a compression. These compressions are sufficiently strong to 
produce a pair of weak shocks in the vicinity of the black hole. These shocks are roughly aligned with 
the line-of-nodes between the black-hole symmetry plane and disk midplane, with one shock above 
the line-of-nodes on one side of the black hole and the other below on the opposite side. These shocks 
enhance angular momentum transport and energy dissipation near the hole, forcing some material to 
plunge toward the black hole from well outside the innermost stable circular orbit. A new, extended 
simulation, which was evolved for more than a full disk precession period, allows us to confirm that 
these shocks and the previously identified "plunging streams" precess with the disk in such a way as 
to remain aligned relative to the line-of-nodes, as expected based on our physical understanding of 
these phenomena. Such a precessing structure would likely present a strong quasi-periodic signal. 
Subject headings: accretion, accretion disks — black hole physics — galaxies: active — MHD 
relativity — X-rays: stars 



1. INTRODUCTION 

The purpose of this paper is to present a more detailed 
analysis of the nume rical simulation results described in 



Fragile et al.| ( 2007b I (hereafter Paper I) . That paper de- 
scribed a global numerical simulation of an accretion disk 
subject to the magneto-rotational instability (MRI) that 
was misaligned (tilted) with respect to the rotation axis 
of a modestly fast rotating (a/M = 0.9) Kerr black hole. 
It was the first such numerical simulation of a tilted disk 
to fully incorporate the effects of the black hole space- 
time and not rely on an ad hoc prescription of angular 
momentum transport. 

A clear identification of a tilted black-hole accretion 
disk in Nature has yet to be made. Nevertheless, there 
is reason to believe they may be quite common. In any 
accretion disk system, the orientation of the disk is set 
by the net angular momentum of the gas reservoir on 
large scales. The orientation of the black hole, on the 
other hand, can either be set by its formation or its evo- 
lution. For stellar-mass black holes in low-mass binary 
systems, only the formation is likely to matter (any evo- 
lution subsequent to the formation is unlikely to change 
the orientation of the black hole significantly) . Since the 
formation of the black hole is largely independent of the 
angular momentum of the gas reservoir, alignment should 
generally not be expected in this case. For high mass bi- 
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naries and active galactic nuclei (AGN), on the other 
hand, a large episode of misalig ned accretion could re- 



orien t the spin of the black hole (Natarajan & Armitage 
19991. Therefore a more detailed understanding of the 
evolutionary history of the system, including mergers in 
the case of AGN, would be needed to know if a tilted 
configuration is expected. 

A tilted disk is subject to differential warping due to 
the effect of Lense-Thirring precession. In di sks with 
a large "viscous" stress (parametrized by the Shakura 



& Sunyaev ( 1973 ) a parameter) the warping results in 



the Bardccn-P etterson configuratio n (Bardeen & Petter- 
1975 Kumar & Pringlc 19851, characterized by an 



son 



alignment of the disk with the equatorial plane of the 
black hole inside some characteristic warp radius. The 
Bardeen-Petterson effect has been invoked by a num- 
ber of authors to explai n peculiar observations such as 
misaligned jets in AGN ( |Kondratko et aL 2005; C apronij 
et al.|2006| [2007 1) and X-ray binaries fragile et al.|2001| 
Maccarone|[^no^| . 

However, we did not see evidence for the Bardeen- 
Petterson effect in our simulation in Paper I. This was 
not surprising since our simulation was carried out in the 
low stress regime, with a < H/r, where H is the half- 
height of the disk. In t his limit warps produced in the 
disk propagate as waves ( |Papaloizou fc Lin|1995 1 , rather 
than diffusively as in the Bardeen- Petterson case. In- 
stead of a smooth transition between an untilted disk 
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at small radii and tilted disk at large radii, as with the 
Bardeen-Petterson effect, we found the tilt to be a long- 
wavelength oscillatory function of radius. 

The tilted simulation in Paper I also showed other dra- 
matic differences from comparable simulations of untiltcd 
disks. Accretion onto the hole occurred predominantly 
through two opposing plunging streams that started from 
high latitudes with respect to both the black-hole and 
disk midplancs. These plunging streams also started 
from a larger radius than the innermost stable circular 
orbit (1SCO), which is often assumed to represent the in- 
ner edge of untilted disks. We interpreted this as being 
due to the fact that the tilted disk encounters a general- 
ized 1SCO sur face at a larger cyl indrical radius than an 
untilted disk (Fragile et al. 2007a I. In this regard the 
tilted black hole effectively acts like an untilted black 
hole of lesser spin. 

Because of the fast sound-crossing time in the disk, 
the torque of the black hole acted globally rather than 
differentially. Instead of strongly warping the disk, the 
torque caused the disk to experience global (solid-body) 
precession. The precession had a frequency of ^ prcc — 
3(Mq/M) Hz, a value consistent with many observed 
low-frequency quasi-periodic oscillations (QPOs). How- 
ever, this value is strongly dependent on the size of the 



disk (Vp rC c oc r , where r a is the outer radius), so this 
frequency may be expected to be vary over a large range. 
In the limit of a very large disk or in a case where the disk 
is strongly coupled to a gas reservoir at large radii, this 
precession frequency may be expected to drop to zero. 

In the present work we expound on some of the features 
of the tiltcd-disk simulation which were not described in 
Paper I. This paper is organized as follows: In Sj2] we 
briefly review the details of the simulations presented in 
Paper I, in particular focusing on models 90h (untiltcd) 
and 915h (tilted), which are identical other than the ini- 
tial tilt of the black hole relative to the disk (/3 = and 
15°, respectively), and model 915m, which is a lower res- 
olution version of 915h. In ^3] we describe the large-scale 
epicyclic motion that appears in the tilted simulation but 
is absent in the comparable untilted simulation. This 
epicyclic motion is driven by radial pressure gradients 
associated with the stationary bending wave created by 
the warping act ion of the gravitomagnetic torque of the 
black hole (e.g. |Nelson fc Papaloizo"u*||1999[ ). Next, in ||4j 
we describe a pair of standing shocks that again have no 
counterparts in untilted simulations. The shocks form 
roughly along the line-of-nodes between the disk mid- 
plane and black-hole symmetry plane, one shock on each 
side of the black hole. We verify that this is not a chance 
alignment by following the long-term evolution of model 
915m and confirming that the shocks remain aligned with 
the line-of-nodes as it precesses with the disk. Finally, 
in ^j5] we consider some of the implications of our dis- 
coveries, particularly in the context of rclativistic iron 
lines. 

2. NUMERICAL METHODS 

The simulations in Paper I were carried out us- 
ing the Cosmos+ + astrophysical magn etohydrodynam- 
ics (MHD) code ( | Anninos et al.||2005| . Cosmos++ in- 
cludes several schemes for solving the GRMHD equa- 
tions; in Paper I, the artificial viscosity formulation was 
used. The magnetic fields were evolved in an advection- 



split form, while using a hyperbolic divergence cleanser 
to maintain an approximately divergence-free magnetic 
field. The GRMHD equations were evolved in a "tilted" 
Kerr-Schild polar coordinate system (t, r, t9, (p). This co- 
ordinate system is related to the usual (untilted) Kerr- 
Schild coordinates (t, r, 9, (f>) through a simple rotation 
about the yy-axis b y an angle (3q, as described in |Fragile] 
& Anninos (|2005 see also Fragile & Anninos 2007j^ 

The simulations were carried out on a spherical po- 
lar mesh with nested resolution layers. The base grid 
contained 32 3 mesh zones and covered the full Air stera- 
dians. Varying levels of refinement were added on top of 
the base layer; each refinement level doubling the resolu- 
tion relative to the previous layer. The main simulations, 
models 90h and 915h from Paper I, had two levels of re- 
finement, thus achieving peak resolutions equivalent to 
a 128 3 simulation. To demonstrate convergence, in Pa- 
per I we also presented results at higher and lower res- 
olutions. In this paper we also discuss results of model 
915m, which used a single level of refinement and had a 
peak resolution equivalent to a 64 3 simulation, but was 
run for significantly longer. 

In the radial direction a logarithmic coordinate of the 
form r\ = 1.0 + ln(r/reH) was used, where tbh = 1-43vg 
is the black-hole horizon radius and tq — GAI/c 2 is 
the gravitational radius. The spatial resolution near 
the black-hole horizon was Ar w 0.05r(j; near the ini- 
tial pressure maximum of the torus, the resolution was 
Ar w 0.57' q. In the angular direction, in addition to the 
nested grids, a concentrated latitude coordinate x-i of the 
form d — X2 + ^(1 — h) sin(2a;2) was used with h = 0.5, 
which concentrates resolution toward the midplane of the 
disk. As a result r cen terAi9 = O.Src near the midplane 
while it is a factor of ~ 3 larger for the fully refined zones 
near the pole. This grid is shown in Figure 1 of Paper I. 

The simulations started from the analytic solution 
for an axisymmetri c torus around a rotating black hole 
(Chakrabarti 1985). The initial torus was identical to 



model KDP of De Villiers et al. (2003), which is the rel- 



ativistic analog of model G'1'4 of |Hawley ( |2000[ ). As in 
model KDP, the spin of the black hole was a/M = 0.9; 
the inner radius of the torus was r- ln — VSro; the radius of 
the initial pressure maximum of the torus r cen ter = 25rg; 
and the power-law exponent used in defining the initial 
specific angular momentum distribution was q = 1.68. 
An adiabatic equation of state was assumed, with T = 
5/3. The torus was seeded with a weak dipole magnetic 
field in the form of poloidal loops along the isobaric con- 
tours within the torus. The field was normalized such 
that initially f3 mag — P/Pb > /?mag,o = 10 throughout 
the torus. For the tilted simulations (915h and 915m) 
the black hole was inclined by an angle Pq = 15° relative 
to the disk (and the grid) through a transformation of 
the Kerr metric. From this starting point, simulations 
90h and 915h were allowed to evolve for a time equiva- 
lent to 10 orbits at the initial pressure maximum, r con t cr , 
corresponding to hundreds of orbits at the ISCO. Model 
915m was evolved for 100 orbital times at the initial pres- 
sure maximum. Table 1 summarizes the parameters of 
all three simulations. 

3. EPICYCLIC MOTION 

As noted above, Lense-Thirring precession causes dif- 
ferential warping in tilted black-hole accretion disks. In 
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TABLE 1 
Simulation Parameters 



Simulation 


Tilt 
Angle 


a/M 


Equivalent 

Peak 
Resolution 


^center a 


Duration 15 


90h 





0.9 


128 3 


25 r G 


10 


915m 


15° 


0.9 


64 3 


25 r G 


100 


915h 


15° 


0.9 


128 3 


25 r G 


10 



a Radius of initial pressure maximum. In units of r. olr b, the geodesic 
orbital period at the initial pressure maximum r con t cr . 



the thick-disk regime, appropriate for the simulations de- 
scribed in Paper I, warp disturbances are expected to 
propagate in a wave-like, rather than diffusive, manner. 
For nearly Keplerian disks, such as those resulting from 
MRI turbulent simulations, the resulting bending wave 
is expect ed to produce large horizon t al motion within 



the d isk (Nelson & Papaloizou 1999 Torkelsson et al. 
2000 1 . The radial and azimuthal components of this mo- 
tion are odd functions of z, leading to large vertical shear 
(dV r /dz, dV*/dz). 

Figures [T] and [2] compare the radial and azimuthal mo- 
tions of the gas for our tilted and untilted simulation 
(915h and 90h). Clearly the fluid velocity in the tilted 
disk is no longer d ominated by turbulent motio n as it is 
for untilted disks (De Villiers & Hawley 20031, but has 
an ordered sense about it. For instance, in Figure [Tja. we 
see that in the right-hand side of the image the gas in 
the upper layers of the disk is moving radially outward, 
whereas the gas in the lower layers is moving radially in- 
ward. The sense of motion is reversed on the left-hand 
side of the image. Similarly, on the right-hand side of 
Figure ^jp, the gas in the upper half of the disk is mov- 
ing slower than the bulk angular velocity, whereas the 
gas in the lower half is moving faster than this average. 
Again, the sense of the flow is reversed on the left-hand 
side of the image. This pattern of motion is reminiscent 
of epicyclic motion, with the upper and lower halves of 
the disk executing epicycles that are 180° out of phase 
with one another. Notice that no such organized motion 
is apparent in our untilted simulation (Figures tth and 

We find that the pattern of the epicyclic motion is tied 
to the orientation of the disk relative to the black hole. 
As our simulated tilted disk 915m precesses, the velocity 
pattern of the epicyclic motion changes as represented 
in Figures [3]and|4j whe re we show the velocity patterns 
at t = 10 and t = 50i O rb, approximately 1/2 precession 
period apart. Note that the sense of the epicyclic motion 
is reversed between the two different evolution times. 

Figure [5] shows the epicyclic motion caused by the 
bending wave from a different perspective. In the figure, 
a 4 x 4 lattice of streamlines starts in the rrz-planc. The 
lattice is centered above and below the symmetry plane 
of the black hole and the viewer is looking almost down 
the black-hole spin axis (the black-hole spin axis is tilted 
10° away from the viewer's line-of-sight to give a bet- 
ter perspective). As expected from the previous figures 
of simulation 915h at t = 10i or b, the streamlines that 
begin above the symmetry plane of the disk are initially 
moving radially outward, whereas those that begin below 



are moving radially inward. After approximately one- 
quarter of an orbit, the upper streamlines have reached 
their apocenter and the radial motion changes direction. 
This is also where material encounters one of the stand- 
ing shocks, described in the next section, which explains 
the sudden change in direction of the streamlines. Much 
of this material then begins plunging toward the hole, 
ultimately accreting within a cou ple orbits. Thus, as 
noted in (Fragile fc Anninos ( 2005 1 , high-latitude mate- 
rial is preferentially being drained from the disk. 

Figure [5] also reveals that the eccentricity e of the par- 
ticle streamlines depends on the height of the streamline 
above the disk midplane (the top streamlines are more 
eccentric that those one row below ). This is consistent 
with the expectation that e cx £ dlvanov & Illarionov 



19971, where £ is the vertical distance measured perpen 



dicular to the midplane of the disk (£ is the z-componcnt 
of a cylindrical coordinate system (i?, ip, £) aligned with 
each concentric ring of the disk, what is sometimes called 
the twisting coordinate system). 



4. STANDING SHOCKS 



The full expression for the eccentricity of the orbit of 
each fluid element is ( Ivanov fc Illarionov||1997 ) 



6M 



6M OR 



(1) 



where (3 and 7 are the tilt and twist of the disk, re- 
spectively, defined for each concentric ring. In Figure [6] 
we plot (3{r), 7(r), and r^/(r) (for computational conve- 
nience we have replaced the cylindrical twisting coordi- 
nate R with the spherical-polar radius r) . It is important 
to note the increase in r\& with decreasing r. Such be- 
havior, with orbits getting more eccentric closer to the 
black hole, results in a concentration of fluid clement tra- 
jectories near their respect ive apocenters (see Fig. 5a of 
Ivanov & Illarionov 19971. Notice in Figure [| for 111- 



stance, that each 01 the top four streamlines makes its 
closest approach to the next streamline out when it is at 
its own apocenter. 

The crowding of fluid clement trajectories near their 
apocenters produces a local density enhancement. This 
will be most pronounced away from the disk midplane 
because of the dependence of e on £. Because the fluid el- 
ements are traveling supersonically, this local density en- 
hancement produces a weak shock in the flow. In Figure 

tve identify the standing shocks by the sudden change in 
e magnitude and direction of the velocity vectors along 
a roughly linear feature associated with the leading edge 
of the each plunging stream. 

Another way to identify the standing shocks is by 
plotting the magnitude of the vorticity, indicated by 
Jcurl V| = |V x V|, which increases at a shock. Figure 
M shows an isosurface of |curl V| overlaid on a density 
isosurface. Very similar results are obtained if the gradi- 
ent of the entropy, another indicator of the presence of 
a shock, is plotted instead. Notice that one of the shock 
surfaces lies mostly above the chosen isodensity surface, 
while the other lies below it. This is consistent with 
the linear dependence of e on £. We don't expect the 
shocks to extend into the disk midplane where epicyclic 
motion ceases. The location of the shocks in azimuth is 
also consistent with the apocenters of the fluid elements 
trajectories, as inferred from Figures [T] — [5j 
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Fig. 1. — Meridional plots (ip = 0) through the final dumps (t = 10t or b) of simulations 915h (left) and 90h (right) showing a pseudocolor 
representation of V for outflowing (V r > 0) and inflowing (V r < 0) gas as hot and cold colors, respectively. The velocity scale is normalized 
to the speed of light c. The plots are overlaid with isocontours of density at p = 0.4. 0.04, and 0.004p maX: o- Material with p < 0.004p max o 
is excluded from the figure. The figure is oriented in the sense of the grid, such that the black hole is tilted 15° to the left in the left panel. 
By this time in simulation 915h the disk has precessed such that the angular momentum axis of the disk is no longer in the plane of this 
figure. At this point Jdisk nas acquired a positive {/-component and a negative x-component. 




Fig. 2. — Meridional plots (ip = 90°) through the final dumps (t = 10t or b) of simulations 915h (left) and 90h (right) showing a 
pseudocolor representation of (V^ — fi)/fi for superorbital (V^ > fi) and suborbital (V^ < f2) gas as hot and cold colors, respectively, 
where Q = (M/r 3 ) 1 / 2 /[l + a(M/r 3 ) 1//2 ] is the particle orbital angular frequency. The plots are overlaid with isocontours of density at 
p = 0.4, 0.04, and 0.004p max o- Material with p < 0.004p maXi o is excluded from the figure. The figure is oriented in the sense of the grid, 
such that the black hole is tilted 15° away from the viewer's vertical in the left panel. 
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Fig. 3. — Meridional plots (<p = 0) at times t = 10 (left) and 50i or b (right) (representing approximately 1/2 precession period difference) 
from simulation 915m showing a pseudocolor representation of V T for outflowing (V r > 0) and inflowing (V r < 0) gas as hot and cold 
colors, respectively. The color scale and contours are normalized the same as in Fig. [I] The fi gure is oriented in the sense of the grid, such 
that the black hole is tilted 15° to the left. Note that the sense of the radial motion as seen from this fixed viewing direction is reversed 
between the two frames. This is consistent with the epicyclic motion tracking the precession of the disk. 




Fig. 4. — Meridional plots (ip = 90°) at times t = 10 (left) and 50t or b (right) (representing approximately 1/2 precession period difference) 
from simulation 915m showing a pseudocolor representation of (V^ — C) /O for superorbital (V^ > Q) and suborbital (V^ < f2) gas as hot 
and cold colors, respectively. The color scale and contours are normalized the same as in Fig. [2] The figure is oriented in the sense of the 
grid, such that the black hole is tilted 15° away from the viewer's vertical. Note that the sense of the residual azimuthal motion as seen 
from this fixed viewing direction is reversed between the two frames. This is consistent with the epicyclic motion tracking the precession 
of the disk. 




Fig. 5. — Isosurface plot of density (semitransparent blue) plus 
selected streamlines for simulation 915h. The isosurface at p = 
O-lpmax.o i s chosen to highlight the plunging streams. There are 
16 streamlines that begin in a 4 X 4 lattice in the irz-plane. The 
lattice is centered about the symmetry plane of the black hole and 
the figure is oriented looking almost directly down the spin axis 
of the hole (red arrow). The spin axis is actually tilted 10° away 
from the viewer's line-of-sight to offer a better perspective. 




Fig. 7. — Isosurface plot of density (semitransparent blue) plus 
selected velocity vectors for simulation 915h at t = 10t or b- The 
density isosurface is the same as in Fig. pj] The figure is oriented 
looking directly down the spin axis of the hole. Velocity vectors are 
only included for gas within ±0.5rQ of the initial disk midplane to 
prevent overcrowding the image. The sudden change in direction 
and magnitude of the velocity vectors (highlighted by the black 
lines) is indicative of a shock. 
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Fig. 6. — Plot of the tilt (/3)t, twist (7)4, and radial dependence 
of the orbital eccentricity rSP as a function of radius through the 
disk. The data for this plot has been time-averaged from t = 9t or b 
to 10t or b- The initial tilt and twist were 0.2618 and 0, respectively. 



Based on our physical understanding of the plunging 
streams identified in Paper I and the standing shocks 
identified in this paper, we expect both features to main- 
tain a constant orientation vis a vis the line-of-nodes be- 
tween the disk midplane and black-hole symmetry plane. 
To confirm this we have evolved our "medium" resolution 
simulation (model 915m of Paper I) for 100i or b, which is 
more than a full precession period. Figure [9] shows the 
plunging streams and standing shocks at 0.1, 0.35, 0.6, 
and 1.1 precession periods. The orientations of both fea- 
tures closely track the precession of the disk and maintain 
roughly constant alignments relative to the line-of-nodes. 




4.1. Shock Effects 



Fig. 8. — Isosurface plots of density (semitransparent blue) and 
[curl V| (red) for simulation 915h at t = 10t or b- |curl V| is a good 
tracer of the location of a shock. The plot is restricted to the region 
7r/4 < 1? < 37r/4 to prevent overcrowding of the image from shocks 
associated with the outflowing jets. The density isosurface is the 
same as in Fig. [5] The figure is oriented looking directly down the 
spin axis of the nole. 



The presence of non-axisymmetric standing shocks 
in the accretion disk can have important conse- 
quences, principally including enhanced angular momen- 
tum transport and dissipation. To illustrate the for- 
mer, in Figure [10] we plot angular momentum resid- 
ual weighted by the circular orbit value, i.e. we plot 
{{{(.) a) t ~ ^chO/^cir &s a function of radius, where £ c - lr is 
the angular momentum of circular orbits with inclina- 
tions of 15° and 0° for simulations 915h and 90h, respec- 
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Fig. 9. — Same as Fig. ^1 except at times t = 10 (panel a), 32 
(panel b), 50 (panel c), and 1004 or b (panel d) (representing approx- 
imately 0.1, 0.35, 0.6, and 1.1 precession period) from simulation 
915m. Also, the isodensity surfaces are at p = 0.04, 0.01, 0.02, and 
0.02p max o, respectively, instead of 0.1p maXj rj as in previous figures. 
The arrows are included to give a visual reference of the precession 
of the line-of-nodes. By noting which shock is above the black-hole 
symmetry plane and which is below we see that: the shocks in 
panels a and b, a and c, and a and d are 90°, 180°, and 360° out of 
phase, respectively. It is apparent that both the plunging stream 
and shock precess with the disk. 

tively, calculated from the following expression (equation 
26 of Paper I) 



ATi + A(Mr) 1 /2iV 2 1 / 2 cos j 



D 



(2) 



where 



iVi 
N 2 



-a 



Mr (3r 2 + a 2 - 4Mr) cos 2 i , (3) 

(4) 



D = a 2 (2r 2 + a 2 



and 



r 4 + a 2 sin 2 i (a 2 + 2r 2 - AMr 

2 ■ . 4 



3Mr) sin 2 i+r 4 +4M 2 r 2 -4r 3 M-Mra 2 

(5) 



A = r 2 - 2Mr + a 2 ; (6) 

and (£)a — (p^)a/(p)a is the density- weighted shell av- 
erage of the specific angular momentum. Shell averaged 
quantities are computed as: 



(SUM) 



2tt r$ 



01 



(7) 



where A = 2 y/—gd'ddip is the surface area of the 

shell. The data have also been time-averaged over the 
interval, 7t rb == tmin < t < t max = 10t or b, where time 
averages are defined as 



(fi)t 



Qdt 



(8) 



The sharp down-turn of the specific angular momentum 
inside r < lOr^ in simulation 915h is indicative of extra 
angular momentum being removed from the flow. This 
suggests that the standing shock plays a significant role 
in the transport of angular momentum in the tilted disk. 
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Fig. 10. — Plot of the weighted residual specific angular momen- 
tum (((£} a) t — ( c ir) / iclr as a function of radius for simulations 915h 
(solid line) and 90h (dashed line). The data has been time-averaged 
over the interval t = 7 to 10t or b- The residuals are calculated from 
the specific angular momenta of circular orbits with inclinations of 
15° and 0°, respectively. 
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Fig. 11. — Plot of the density- weighted shell average of the fluid 
entropy as a function of radius for simulations 915h (solid line) 
and 90h (dashed line). The data has been time-averaged over the 
interval t = 7 to 10t or b- 

The shock also enhances the energy dissipation in the 
disk. This is shown in Figure where we plot the 
density-weighted shell averages of the fluid entropy in 
the tilted and untilted disks. Actually, since we are only 
interested in relative changes in entropy, for simplicity 
we plot S — ln(P/p r ). The significant enhancement in 
entropy generation inside r < lOr^ for the tilted disk 
indicates the extra dissipation provided by the shock. 
This up-turn nicely coincides with the down-turn in Fig- 
ure 10 further supporting the association of both effects 
withthe shock. 

The extra energy dissipation and angular momentum 
transport significantly affect the inner regions of the disk. 
In Figure 1121 we see that the radial plunging region, de- 

T 

fined by the sharp upturn in V , begins at a consider- 
ably larger radius for the tilted simulation (915h) than 
for the untilted (90h). This c an have important impli- 
cations for disk observations ( |Krolik fc Hawley 2002| ) 
since most of the radiated energy from a disk comes 
from just outside the plunging radius. Of special con- 
cern is the common use of the inner edge of the disk as a 
direct indicator of the spin of the black hole. Clearly 
tilt must be taken into account if one wants to re- 
late the plunging radius of the disk to the spin of the 
black hole for tilted disks. Figure 12 also compares 
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Fig. 12. — Plot of characteristic velocities as functions of radius within the tilted (left panel) and untilted (right panel) disks. All data 
have been time-averaged over the interval t = 7 to 10t or b. Of particular importance is the sharp upturn in V r , indicating the start of the 
plunging region, which happens at a larger radius for the tilted simulation. 



other characteristic velocities associated with the disks. 
All of the velocities are density weighted shell averages 
V = ((pV) a I '(p) A)t- The local sound speed is recovered 
from cl = T(r - l)P/[(r - l)p + TP]. The Alfven speed 
is 
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We approximate the turbulent velocity as Vturb 
where 

\u r U v\\B\\ 2 - B r B ip \ 
A-kP 



(10) 



is the dimensionless stress parameter. The most notable 
difference is in the radial inflow velocity V , although 
there is also some indication that the turbulent stress 
rises faster in the interior of the tilted disk. 

4.2. Post-Shock Magnetic Field 

We found previously that magnetic fields remain sub- 
thermal everywhere in the disk throughout the simula- 
tion (see Figure 11 of Paper I). They, therefore, must not 
play a significant role in the dynamics of the epicyclic mo- 
tion, plunging streams, or standing shocks. This state- 
ment is further supported by the fact that very similar 
features were seen in earlier hydrod ynamic simulations 
of til ted black-hole accretion disks ( [Fragile fc Anninos 
2005|) . Despite the minor role of the magnetic fields in 



the dynamics, the presence of a standing shock can en- 
hance the st rength of the magnetic field. From |Fragile| 
et al. (20051 we note that the post-shock value of /3 mag 
in the Newtonian limit can depend sensitively on the 
strength of the shock 

/3 mag ,p S = 2r(r - i) 2 /(r + i) 3 M 2 /3 magii (11) 

where /3 ma g,i is the pre-shock value and we assume the 
field is oriented perpendicular to the shock normal. For 
r = 5/3, as in our work, this gives 

Anag,ps = 0.078X 2 /3 magjl . (12) 
Therefore, for M < 3.6, the post-shock / 3 ma g will be 
lower than the pre-shock value. In Figure |13| we show 
that there is, indeed, a thin layer of magnetically domi- 
nated pla sma just behind the shock. To get /3 ma g,ps = 0.1 
as Figure 13 starting from /3 magi i = 1 (a reasonable guess 
for the highlatitude material), we need M. — 1.13, which 
is consistent with what we sec in the simulation. Again, 
these do not appear to be tremendously strong shocks. 



5. DISCUSSION 

In this paper we have compared two simulations that 
are identical to one another in every respect except one: 
the initial tilt between the black-hole and disk angular 
momenta. Despite the apparent similarity, we observed 
remarkable differences in the evolution of these two sim- 
ulations. 

The primary new feature that we describe in the main 
disk body is a strong epicyclic driving attributable to the 
gravitomagnetic torque of the misaligned (tilted) black 
hole. The induced motion of the gas is coherent over 
the scale of the entire disk. An interesting point about 
this epicyclic motion that has not been made before is 
that it could be detectable, for instance in the profile of 
rclativistically broadened iron lines. We have previously 
pointed out the importance o f the iron lines in directly 
probing tilted accretion disks ( Fragile et al.|2005| , but at 
the time we were not in a position to recognize the im- 
portance of the epicyclic motion. From Figure [2] we get 
that the velocities associated with the epicyclic motion 




Fig. 13. — Isosurface plot of density (semitransparent blue) plus 
an isosurface of /3 ma g = P/Pb = 0.1 (orange), indicating a thin 
layer of magnetically dominated plasma behind the shock. The 
density isosurface is the same as in Fig. The figure is oriented 
looking directly down the spin axis of the hole. [5] 
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represent a significant fraction (< 40%) of the orbital ve- 
locity of the gas. The corresponding shift in a reflection 
feature such as the iron line should be of a similar mag- 
nitude. The interesting thing is to note that an observer 
viewing the disk from the vantage point of Figure [2]a and 
seeing the "top" of the disk would see a smaller than 
expected red shift (the gas going away is not moving as 
fast as expected) and a larger than expected blue shift 
(the approaching gas is moving faster than expected). 
An observer viewing the same disk from the same van- 
tage point but seeing the "bottom" of the disk would see 
exactly the opposite shift. Therefore, depending on the 
viewing angle, the entire line profile of a misaligned disk 
could be shifted toward the blue or the red relative to 
an aligned disk. On the redshifted side this effect might 
be confused with gravitational redshifting, making it ap- 
pear that the line is coming from deeper in the potential 
well than is actually the case. 

If the disk actually precesses then there is a clear and 
simple way to disentangle this effect because the epicyclic 
motion is phased with the orientation of the disk, as dis- 
cussed in f|3] above. This means that the line shift we are 
describing would reverse itself twice per precession pe- 
riod, appearing for half a precession period as an overall 
blueshift and for the other half as an overall redshift. If 
the integration time of the detector is longer than the 
precession period, the net effect would generally be to 
smear or broaden the line. This may provide an alter- 
native explanation for an exceptional ly broad line profile 
such as GX 339-4 ( |Miller et aL"|2004| ), without requiring 
a high black-hole spin or small disk radius. On the other 
hand, if the integration time is shorter than the preces- 
sion period, then one would expect the iron line to vary 
in phase with changes in the X-ray flux (i.e. in conjunc- 
tion with a low frequency QPO that would correspond 
to the precession frequency of the disk) . Such behavior 
has been identified in at least one source, GRS 1915+105 
( |Miller fc Hmrla^|[2005l ). 

The second new feature of our tilted disk simulation 
that we described is a pair of standing shocks, roughly 
aligned with the line-of-nodes between the disk and black 
hole symmetry planes. This asymmetric shock provides 
additional angular momentum transport and energy dis- 
sipation in the tilted disk, relative to the untilted one. 
This enhanced dissipation may help compensate for the 
loss of radiative efficiency (relative to an untilted disk) 
due to the plunging region starting further out than the 
equatorial ISCO radius of the spinning black hole. In 
collisionless accretion fl ows, such as tho se thought to be 
relevant for Sgr A* (e.g. Quataert|2003[ ), it is conceivable 



that such shocks could also be sites of particle accelera- 
tion. The precession of these shocks might then result in 
periodic variations of nonthermal radiation. 

In considering the magnetic fields in the inner part of 
the tilted disk, as noted in Paper I, the magnetic fields 
remain largely subthermal. They do not play an impor- 
tant role in the physics of the epicyclic motion, plunging 
streams, or standing shocks, although we do find some 
enhancement of /3 m ag associated with the shocks. 



We finish with a figure ( 14 1 which provides a visual 
summary of our findings. We note the overall consis- 
tency of our interpretation of the results: 1) The ob- 
served epicyclic motion is in agreement with expecta- 
tions for warped disks in the wave-like propagation limit; 
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shock 
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Fig. 14. — Schematic diagram of the inner region of the tilted ac- 
cretion disk, showing the pattern of epicyclic motion, the standing 
shock, and the plunging streams. 

and 2) The shocks are located near the apocenters of the 
epicyclic motion, as one expects when the eccentricity 
increases with decreasing radius. 



We would like to recognize Christopher Lindner and 
Joseph Niehaus for their contributions to this work. We 
also thank Steve Balbus and Gordon Ogilvic for very use- 
ful discussions and comments and the anonymous ref- 
eree for suggested improvements. PCF gratefully ac- 
knowledges the support of a Faculty R&D grant from 
the College of Charleston and a REAP grant from the 
South Carolina Space Grant Consortium. This work was 
supported under the following National Science Founda- 
tion grants and programs: AST 0707624, Partnerships 
for Advanced Computational Infrastructure, Distributed 
Terascale Facility (DTF), and Terascale Extensions: En- 
hancements to the Extensible Terascale Facility. 



10 



REFERENCES 



Anninos, P., Fragile, P. C, & Salmonson, J. D. 2005, ApJ, 635, 723 

Bardeen, J. M., & Petterson, J. A. 1975, ApJ, 195, L65 

Caproni, A., Abraham, Z., Livio, M., & Mosquera Cuesta, H. J. 

2007, MNRAS, 379, 135 
Caproni, A., Abraham, Z., & Mosquera Cuesta, H. J. 2006, ApJ, 

638, 120 

Chakrabarti, S. K. 1985, ApJ, 288, 1 

De Villiers, J., & Hawley, J. F. 2003, ApJ, 592, 1060 

De Villiers, J., Hawley, J. F., & Krolik, J. H. 2003, ApJ, 599, 1238 

Fragile, P. C, & Anninos, P. 2005, ApJ, 623, 347 

Fragile, P. C, Anninos, P., Blaes, O. M., & Salmonson, J. D. 

2007a, in proceedings of the 11th Ma rcel Grossmann Meeting 

on General Relativity (|astro-ph /0701272 ) 
Fragile, P. C, Anninos, P., Gustalson, K., & Murray, S. D. 2005, 

ApJ, 619, 327 

Fragile, P. C, Blaes, O. M., Anninos, P., & Salmonson, J. D. 2007b, 
ApJ, 668, 417 

Fragile, P. C, Mathews, G. J., & Wilson, J. R. 2001, ApJ, 553, 955 
Fragile, P. C, Miller, W. A., k, Vandernoot, E. 2005, ApJ, 635, 157 



Hawley, J. F. 2000, ApJ, 528, 462 

Ivanov, P. B., & Illarionov, A. F. 1997, MNRAS, 285, 394 
Kondratko, P. T., Greenhill, L. J., & Moran, J. M. 2005, ApJ, 618, 
618 

Krolik, J. H., & Hawley, J. F. 2002, ApJ, 573, 754 
Kumar, S., & Pringle, J. E. 1985, MNRAS, 213, 435 
Maccarone, T. J. 2002, MNRAS, 336, 1371 
Miller, J. M., & Homan, J. 2005, ApJ, 618, L107 
Miller, J. M., et al. 2004, ApJ, 601, 450 
Natarajan, P., & Armitage, P. J. 1999, MNRAS, 309, 961 
Nelson, R. P., & Papaloizou, J. C. B. 1999, MNRAS, 309, 929 
Papaloizou, J. C. B., & Lin, D. N. C. 1995, ApJ, 438, 841 
Quataert, E. 2003, Astronomische Nachrichten Supplement, 324, 
435 

Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337 
Torkelsson, U., Ogilvie, G. I., Brandenburg, A., Pringle, J. E., 
Nordlund, A., & Stein, R. F. 2000, MNRAS, 318, 47 



